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Abstract. We attempt to constrain the luminosity function of Gamma Ray Bursts (GRBs) from the observed 
number count -flux relation and the afterglow redshift data. We assume three classes of luminosity functions for 
our analysis: (a) Log-normal distribution, (b) Schechter distribution, and (c) Scale-free distribution. We assume 
several models of the evolution of the GRB population for each luminosity function. Our analysis shows that: 
(a) log-normal is the only luminosity function that is compatible with both the observations. This result is 
independent of the GRB evolution model, (b) for log-normal function, the average photon luminosity Lo and the 
width of the luminosity function a that are compatible with both the observations fall in the range: 10 55 sec -1 
Lo Si 10 56 sec -1 and 2 a £ 3, (c) the agreement of observations with other luminosity functions requires the 
GRB population to evolve more strongly than the evolution of the star-formation rate of the universe. 
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1. Introduction 

Recent afterglow observations of GRBs have given im- 
portant clues about the cosmological nature of GRBs 
and their environment and luminosities (for details see 
Kulkarni et al. 2000 and references therein). The afterglow 
database (see Greiner's home page; Greiner (2000)) is suf- 
ficiently large that one could now think of using these data 
for other studies in astronomy. For instance, one could ad- 
dress the issue concerning the luminosity distribution of 
GRBs using the redshift measurements available in the 
database. The luminosity distribution of GRBs was hith- 
erto obtainable only from the number count-flux relation- 
ship (Piran 1992, 1999; Mao & Paczynski 1992; Woods & 
Loeb 1995, Lubin & Wijers 1993; Ulmer & Wijers 1995; 
Ulmer et al. 1995; Cohen & Piran 1995, Band, Hartman, 
& Schaefer 1999, Hakkila et al. 1996, Horack & Hakkila 
1997). The redshift measurement of several GRB after- 
glows has allowed direct determination of GRB luminosi- 
ties. The distribution of these luminosities can indepen- 
dently be used to get information about the luminosity 
function. 

In this paper we constrain the luminosity function of 
GRBs using two approaches: (a) using the redshift data 
for a sample of 16 GRBs from afterglow measurements; 
(b) using the number-count v/s flux (i.e, Af—F ) rela- 
tion for GRBs in the current Burst and Transient Source 
Experiment (BATSE) catalog. In doing so, we consider 



various luminosity functions (viz. Schechter, scale-free and 
log-normal) each with 'no- evolution' and some of the evo- 
lutionary models consistent with star-formation history of 
the universe. Our aim is to identify the class of luminosity 
functions which are consistent with both the Af—F rela- 
tion for BATSE GRBs and the luminosity distribution of 
the sample of GRBs whose redshifts are available from the 
afterglow observations. 

In §2 we briefly review the Af—F method and de- 
scribe models of luminosity function and GRB number 
density evolution. In §3 we discuss the afterglow redshift 
data and the Likelihood method used for extracting in- 
formation about the luminosity function parameters from 
this data. The results are presented and discussed in §4. 
Throughout this paper we assume a cosmological model 
with Q m = 0.3 ,n A = 0.7 ,H Q = eSkmsec^Mpc -1 ; this 
model is favoured by recent observations (Perlmutter et 
al. 1998, de Bernardis et al. 2000, Freedman et al. 2000). 



2. The M(>F)-F of BATSE GRBs 

For calculating the number of bursts exceeding a given 
flux F, Af(>F), from BATSE sources we first need to de- 
fine 'flux' of a GRB. We use the BATSE peak photon flux 
(Paczynski 1995 ) averaged over the trigger time of 1.024 
seconds as being representative of GRB flux. BATSE re- 
ports peak photon flux, F, integrated over an energy range 
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from 50 kev to 300 kev . This is related to the photon lu- 
minosity L (in photons sec^ 1 ) as (see Appendix A for a 
derivation) : 



L(l + z) 2 
4nDl 



(1) 



Here F is in units cm _2 sec _1 . a is the spectral index of 
intrinsic photon luminosity. We take a = 2 as it provides 
a reasonable ht to the burst spectra (for more details and 
caveats see Band et al. 1993. Note that the spectral index 
a here is the high-energy spectral index /3 as defined in 
Band et al. 1993). The luminosity distance Dl = r(l + z), 
r is the coordinate distance to the GRB at a redshift z. 
For flat cosmologies i.e. f^Totai = 1, it is given by: 



r = H~ 



dz' 



o (fi m (l + z'f + Oa) 1 / 2 



(2) 



Here the Hubble length Hq 



9.1 x KF/i" 1 cm. n Ti 



and Oa are the present energy densities (in units of the 
critical density) of the non-relativistic matter and the 
cosmological constant, respectively. Observations of SNIa 
at high redshifts and anisotropics of Cosmic Micro-wave 
Background Radiation (CMBR) are consistent with a flat 
cosmological model with fl\ ~ 0.7 (Perlmutter et al. 1998; 
de Bernardis et al. 2000). 

In current BATSE catalog peak fluxes are available for 
2093 events. We take F > 0.4 cm~ 2 sec -1 for comparing 
theoretical predictions with the observed AT(>F)-F from 
the BATSE catalog (Loredo & Wasserman 1998). This 
requirement leaves ~ 1790 burst for our analysis. These 
GRBs are divided into 170 flux bins for computing the 
number count- flux relation. 

The number count Af(>F) can be expressed as: 



JV(>F) = 4tt 



Here 



dr~ 2 



v- I dL n{L,z)(l + z)- a .(3) 

4ttD? F 



n(L,z) — <f>*(z)(f)(L)dL 



(4) 



is the comoving number density of GRBs in a given lumi- 
nosity range L to L + dL at a given redshift. <j>(L) is the 
luminosity function, defined such that J Q </)(L)dL = 1. 
4>*{z) — </>*(0) x (1 + z) 1 is the total comoving number 
density of GRBs at redshift z. The factor of (1 + z)~ a in 
Eq. (^) gives the 'k-correction' which corrects for the dif- 
ference between the emitted and the observed wavelength. 

The functional form of the luminosity function of the 
GRBs is unknown. In the framework of fireball models, the 
luminosity function of GRBs is expected to be broad with 
luminosity width of nearly 2 orders of magnitude (Kumar 
& Piran 1999). For our study, we assume several different 
luminosity functions, <f>(L): 

Log-Normal distribution function: 



exp(- ( i 2 /2) 

<^ L ) = K 9 

V 27TCH 



-(ln(L/L )) 2 /(2^ 2 )] ±- (5) 



er and Lq are the width and the average luminosity of 
the luminosity function, respectively (for previous use of 
this luminosity function for GRB studies see e.g. Woods 
& Loeb 1995). The log-normal distribution is also repre- 
sentative of the luminosity function of spiral galaxies (for 
details see Bingelli, Sandage & Tamman 1988). 
Schechter distribution: 



-0 



exp (— L/L*) 



1 



(6) 



A is a normalizing constant. The entire galaxy population 
can be roughly represented by this function (for details 
see Bingelli et al. 1988). We restrict ourselves to f3 < 1 in 
this paper. 

Scale-free luminosity function: 



0(L) = A [ — — for L B 



< L< Lr, 



(7) 



A is a normalizing constant. Several authors have used 
this class of luminosity functions for GRB analyses (see 
e.g. Schaefer 2000). 

Af(>F) depends on a number of parameters: (a) cos- 
mological parameters through Dl and r. We fix the cos- 
mological parameters to their most favoured values, (b) 
Zmax, the maximum redshift of GRBs, (c) <f>*{z), which 
gives the redshift dependence of the GRB population, (d) 
the parameters of the luminosity function. Our aim is to 
obtain the most-favoured values of the luminosity function 
parameters. After fixing the cosmological model and the 
GRB spectral index, the most important remaining uncer- 
tainty comes from the evolution of GRB population. One 
obvious choice is the 'no evolution' model. In this model, 
(j>*{z) is independent of z. However afterglow observations 
give circumstantial evidence that GRBs might be associ- 
ated with star-forming regions (Kulkarni et al. 2000 and 
references therein). Such an association might mean that 
GRB population trace the star-formation history of the 
universe. Therefore we also consider several other models, 
which are consistent with the star-formation history of the 
universe: 



Model I: No Evolution model 

Model II: The GRB population evolves as the luminos- 
ity density at 1600 A, which is a good tracer of star- 
formation history (Madau et al. 1998, Madau et al. 
1996). In this case the comoving number density of 
GRBs increases by a factor of ~ 10 from z ~ to 
z ~ 1.5 (Lilly et al. 1996) and then decreases approx- 
imately oc (1 + z)^ 1 at higher redshifts (Figure 3 of 
Madau et al. 1998). 

Model III: To account for the possible dust-extinction in 
the star- forming regions for z > 1.5 (Conolly et al. 



1 for more details see [http : //gammaray .msf c .nasa. gov/batse/datdi| 97), we assume that the star-formation rate grows 
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for z < 1.5 as in the previous model, but remains con- 
stant at the value of z — 1.5 till z ~ 3, and then 
declines as oc (1 + z) 1 at larger redshifts. 

Model IV: We consider an extreme model that, given the 
uncertainty in the evolution of the star-formation rate 
at z > 1.5, is barely consistent with the star- formation 
history of the universe. In this model, the comoving 
number density of GRB population evolves oc (1 + z) 3 5 
for z 1.5 (Lilly et al. 1996). The comoving number 
density remains at the value of z = 1.5 for z < 5 
and then it gradually declines at oc (1 + z)~°' 5 . In 
this model much of star- formation at redshifts 1.5 
occurs in regions highly shrouded by dust which al- 
lows a high star-formation rate to be compatible with 
the star-formation rate observed in the Hubble Deep 
Field UV drop-out galaxies. However, further increase 
in the star-formation rate, which must be accompa- 
nied by a suitable increase in the dust content, might 
be incompatible with the observation of Far Infra-red 
background (Puget et al. 1996, Guiderdoni et al. 1991). 

3. Observed GRB redshifts 

In Table 1 we list the redshifts of GRBs used to calcu- 
late the luminosity distribution of GRBs. The GRB at 
z = 0.008 (GRB 980425) has been excluded from this list, 
as it is probably associated with a supernova (Galama et 
al. 1998) and therefore corresponds to a different pop- 
ulation of GRBs (Kulkarni et al. 2000). Three GRBs 
whose redshifts are not certain (GRB 980326: z ~ 1; GRB 
980329: z < 3.5; GRB 990507: z ~ 0.25) are also excluded. 

Gamma-ray fluxes (photons cm -2 sec -1 ) correspond- 
ing to these 16 GRBs arc listed in Table 1: Eight of these 
are taken from BATSE flux table (integrated over 1024 
milliseconds). The gamma-ray fluxes of remaining GRBs, 
triggered either by BeppoSAX or IPN satellites are in dif- 
ferent energy band and therefore had to be extrapolated 
to energy range of BATSE. The band-pass conversion is 
achieved using photon luminosity spectral index, a = 2 
(see appendix B). We adopt the values from Schaefer 
(2000) for three of the GRBs. 

From the measured redshifts and observed 7— ray 
fluxes, the total energies of GRBs, assuming isotropic 
emission, is in the range 10 51-54 erg. The production of 
such large energies pose serious problems in theoretical 
modelling. One of the proposal to resolve the energy cri- 
sis was to have the emission collimated into jets, which 
would lower the total energies by a factor of few hun- 
dreds (Rhoads 1999; Piran 1999). The GRB afterglow ob- 
servations have shown evidences for beaming in the light 
curves in 5 cases : GRB990123 (Holland 2000; Kulkarni 
1999), GRB990510 (Holland 2000, Harrison 1999), GRB 
991216(Halpern 2000), GRB 000301c (Sagar et al. 2000, 
Berger et al. 2000 ) and GRB 000926 (Price et al. 2000). 
The opening angle 9 of the jet may be calculated from 
Eq. (1) of Sari et al. (1999) knowing the break in light 
curve. The beaming factor in Table 1. is ~ 2/6* 2 and has 
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been used to apply the corrections in luminosities in 5 
cases. 

The probability that burst of a given flux will occur in 
a redshift range z to z + dz is given by: 

p(z)dz = 4>{L) x -j-dz (8) 

4>{L) is given by Eqs. (||), or ((?]), and L for a given 
flux F and redshift is determined by Eq. (|l|). 

The likelihood that the observed GRB luminosities 
were drawn from a given luminosity function is: 

16 

C(a k )=l[p(z i) F i ) (9) 

We maximize the likelihood function with respect to the 
parameters of the luminosity function, a*,, e.g. = 
{Lo,a} for the Log-normal luminosity function. It should 
be noted that parameters estimated using this method do 
not depend on the number density <j>, (z) of the GRBs. 

4. Results 

We show our results for the log-normal luminosity function 
in Figure 1. We plot the allowed region in the Lq-u plane 
from both the observed M{>F)-F from BATSE catalog 
and the afterglow redshift data. The allowed region from 
requiring consistency between the observed Af(>F)-F re- 
lation and the theoretical number counts (Eq. ^) corre- 
sponds to the area in which the K-S probability Pks > 0.01 
(see e.g. Press et al. 1992 for details on the K-S test). 

The best fit value of a and Lq from the afterglow red- 
shift data are: a — 1.7 and Lq = 3x 10 57 sec" 1 without the 
beaming corrections and a = 2 and Lq = 2 x 10 56 sec _1 
with the beaming corrections. We show in Figure 1 the 
region within which the value of the likelihood function is 
> 10~ 4 times the value at the maximum. This correspond 
roughly to 99% confidence level for a two-parameter fit 
(Press et al. 1992). (We do not calculate the joint con- 
fidence levels using the usual Fisher matrix approach be- 
cause the Likelihood function is not a joint Gaussian dis- 
tribution and also because the Likelihood function is very 
broad near the maximum). 

Results for other assumed luminosity functions are 
shown in Figure 2, 3 and 4, using the same criteria for 
showing the allowed regions as given in the preceding para- 
graph. The results are shown for only two representative 
scale-free models. The best-fit values of parameters from 
GRB redshift data are: (3 = 0.6, L* = 3 x 10 59 sec -1 with- 
out beaming correction and (3 = 0.6, I, = 2 x 10 59 sec _1 
with beaming correction (Schechter luminosity function); 
(3 = 0.65, L* = 5 x 10 56 sec _1 without beaming correc- 
tions and (3 = 1.05, L* = 3.8 x 10 56 sec _1 (scale-free 
model with L m - m ~ and L max — 10 3 £»); f3 — 0.65, 
L, = 6 x 10 57 sec _1 without beaming correction and 
(3 = 1.15, L* = 1.2 x 10 58 sec _1 with beaming cor- 
rection (scale-free model with L min ~ 3 x 10~ 2 L* and 
Lmax — 10 2 .L*). 
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As seen in the Figures, it is possible to explain the 
observed Af(>F)-F relation using any of the luminosity 
functions we assume. Further the allowed regions do not 
strongly constrain either the luminosity of the GRB or the 
width of the luminosity function (Loredo & Wasserman 
1998, Schmidt 2000). For any luminosity function the al- 
lowed range of luminosities span atleast a decade depend- 
ing on the evolution of GRB population (Figure 1). It 
can be noticed that stronger evolution in the GRB pop- 
ulation results in a higher average luminosity (Schaefer 
2000). z mal , the maximum redshift of GRBs is not an 
important parameter so long as z max <; 5. It is be- 
cause even the faintest bursts we consider in our analysis 
(F ~ 0.4 cm -2 sec -1 ) come from z ^ 5 for much of the 
allowed luminosity range. 

The GRB redshift data give a large range of photon lu- 
minosities (10 57 sec -1 ^ L ^ 10 59 sec -1 in photons sec -1 ). 
Therefore it is natural to expect that the underlying lu- 
minosity function is broad. It is clearly seen in the figures. 
It should be noted that the GRB luminosities implied by 
the Af(>F)-F relation are not at variance with the ob- 
served GRB luminosities. However, as seen in the figures, 
these two observations imply two different set of luminos- 
ity function parameters in most cases. 

The log-normal is the only luminosity function that is 
compatible with both BATSE number count and the af- 
terglow redshift data. This agreement is independent of 
the GRB evolution model. The agreement between the 
two data sets requires 10 55 sec -1 ^ L ^ 10 56 scc -1 and 
2 £3 a ^ 3. For other luminosity functions the agreement 
between the two observations becomes better as GRB evo- 
lution becomes stronger (Figures 2, 3, and 4). However 
the requirement that the two regions overlap can be ful- 
filled only if the GRBs evolve at a rate much stronger than 
the star-formation rate of the universe. 

It is seen in Figures 1,2,3, and 4 that the beaming cor- 
rection does not make any quantitative difference to our 
results. In almost all the cases the 99% region from the 
afterglow data becomes smaller after the beaming correc- 
tion, but it is difficult to draw any conclusions from it. A 
particularly interesting case is one of the scale-free mod- 
els (Figure 4) in which the beaming correction reduces 
the area of the 99% region quite dramatically. However 
this is owing to the fact that the beaming correction in- 
creases the luminosity width of the observed GRB after- 
glows, and the model in question has a luminosity width 
of 1000 which matches the luminosity width of the GRBs. 
Therefore only of a very small range of L* is compatible 
with observations. 

Selection effects: One possible reason for the dis- 
crepancy/agreement between the number count-flux and 
GRB redshift data results could be selection effects. Most 
of the redshift determinations have been possible owing 
to the detection of X-ray counterpart of the GRB by 
BeppoSAX. This results in two types of selection effects: 
(a) BeppoSAX is not sensitive to burst duration < 1 sec 
(Feroci et al. 1999), (b) the GRBs detected by BeppoSAX 



might be X-ray selected though Feroci et al. (1999) em- 
phasize it is unlikely. 

The burst-duration selection effect might mean the 
GRBs for which the redshifts have been determined belong 
to a different population of GRBs. The burst duration dis- 
tribution for the BATSE bursts is observed to be bi-modal 
(Fishman & Meegan 1995), which strongly suggest that 
long- and short-duration bursts are different populations 
of GRBs. In light of this it is usual to expect that these 
two might have different luminosity functions. To check 
for this effect we make a sub-catalog of GRB bursts with 
duration > 2 sec (measured in T90) . As seen in Figures 1 
to 4, the range of allowed luminosity function parameters 
from the Af(>F)-F relation of this sub-catalog are qual- 
itatively similar to those for the entire catalog, and our 
results arc not sensitive to this selection criterion. 

Another selection effect is Malmquist bias (see e.g. 
Sandage 1993 for a discussion and relevant references): 
If the luminosity function of sources is calculated from a 
magnitude-limited sample then the average luminosity of 
the estimated luminosity function exceeds the true mean 
of the underlying population. Other moments of the lu- 
minosity function (variance, etc.) are also affected. For 
GRBs, unlike other astronomical sources, it is difficult to 
establish the criterion of source selection .i.e. it is not com- 
pletely clear that the source whose luminosities have been 
established are selected randomly from a flux- limited sam- 
ple. For the purposes of this paper we have taken this as- 
sumption to be true. To assess the effect of Malmquist 
bias, we calculate, for the luminosity functions considered 
in this paper, the average luminosity (L) f and its vari- 
ance for a flux-limited sample (see e.g. Sandage 1993). 
The flux limit of the sample is taken to be the mini- 
mum flux for which the luminosity has been determined 
(~ 0.5 cm -2 sec -1 ; Table 1). In addition to flux, (L)f has 
dependence on the parameters of the luminosity function 
and evolutionary properties of the whole population. We 
find that, within the range of relevant luminosity function 
parameters and the evolution of GRB population: 

1. The average luminosity of the flux- limited sample ex- 
ceed the true mean by a factor of two to three. 

2. The variance of the flux- limited sample is within 
40-50% of the true variance. 

This is well within the range of 99% confidence contours 
from the afterglow observations (Figures 1 to 4) . It should 
be pointed out that to get complete information about the 
luminosity function, we need to calculate all the higher 
moments. However, we note that uncertainties in the first 
few moments are far below the range of luminosity func- 
tion parameters allowed by the likelihood analysis of the 
afterglow data (Figures 1 to 4). This suggests that errors 
in the luminosity function parameter are more due to the 
smallness of the sample than the Malmquist bias, and our 
conclusion are not affected by this bias. 

To sum up: of the four luminosity functions consid- 
ered here, log-normal is the only luminosity function that 
is consistent with both the BATSE number count-flux 
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relation and afterglow redshift data, independent of the 
evolution of the GRB population. Notwithstanding the 
selection effects, the other luminosity function can be con- 
sistent with both observations only if the GRBs evolution 
far exceeds the evolution of star-formation rate in the uni- 
verse. 



Appendix A 

Here we give a brief derivation of Eq. (|l|) . Assume that a 
source at a redshift z is emitting with photon luminosity 
L v (photons sec -1 Hz -1 ), with intrinsic spectrum: 



L v = L t ( — 



(10) 



Here vj is some fiducial frequency. The received photon 
flux (photons cm -2 sec -1 Hz -1 ) is: 

F »° = 4^ ("> 

Here r, defined in Eq. (||), is the coordinate distance to 
the source from the observer at the present epoch and Vq 
is the observed frequency. Note that there are no factors 
of (1+z) in the denominator of Eq. (fill). Integrating the 
flux over the observed band-pass (50kev to 300 kev for 
BATSE), and using v/v Q = (1 + z), we get: 



F = F Vo dv = 



Aur 2 



du 



(12) 



(13) 



Defining 

L = Lf I **' 
and using = (1 + z)r, we obtain Eq. ([!]). 

Appendix B 

Throughout this paper we take the 'flux' to mean the pho- 
ton flux averaged over 1.024 sec BATSE trigger, in the 
energy range 50-300 kev. Out of 16 afterglows whose red- 
shifts have been determined, only 8 have known BATSE 
fluxes. In the other cases, the gamma-ray fluxes are ob- 
served in other frequency bands by either BeppoSAX or 
IPN satellites. To be consistent with our definition of flux, 
we extrapolate the observed fluxes from the observed band 
to the BATSE band. We briefly describe the method of 
this extrapolation in this appendix. It should be pointed 
out that the difference of trigger time between BATSE 
and other instruments must also be taken into account for 
extrapolating fluxes. However, we apply only band-pass 
corrections. 

Assume that a source at redshift z is emitting with a 
luminosity L v (erg sec -1 cm -2 ) and spectrum: 



L V = LA — 



(14) 



Vf is some fiducial frequency. The observed flux at the 
frequency = u/(l + z) is: 



F 

47rr 2 (l + z) 



(15) 



The flux integrated over a band-pass between frequencies 
v\ and V2 is given by: 



F(vi, v 2 )= I F yQ dv Q 



(16) 



with the definition F(z) = L f (l + z) -/3-1 /( 47rr2 )- We as- 
sume j3 = 1 in this paper (the spectral index a for the 
photon luminosity is related to (3 as a — (i + 1). This 
gives: 

F(v 1 ,v 2 ) = ^-\og(v 2 /v 1 ). (17) 
v f 

Eq. ( |l7j ) can be used to convert flux observed in any band- 
pass to the BATSE band-pass. 

The photon luminosity (LP) is related to the energy 
luminosity as: Lf, — L v l(hv), h being the Planck's con- 
stant. This relation and the methods described above as 
well as in the previous appendix can be used to get the 
conversion between the band-pass integrated photon flux 
(F) and the energy flux (F '(vi , v 2 )) ■ The relation is (for 
(5 = 1): 



F = F{v 1 ,v 2 ) 



h log(^ 2 M) 



(18) 
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Table 1. Redshifts and luminosities of Afterglow GRBs 



GRB 


BATSE 


Redshift 


Photon flux 


Luminosity 


Beaming 


Ref. for 




tr# 


z 


ph/s/cm 2 


ph/s 


factor 2 * 


photon flux 3 ** 


GRB970228 




0.695 


10.0 


2.34 x 10 bs 




1 


GRB970508 


6225 


0.835 


0.969 


3.56 x 10 s7 




2 


GRB970828 


6350 


0.958 


1.5 


7.74 x 10 57 




1 


GRB971214 


6533 


3.418 


1.955 


2.27 x 10 59 




2 


GRB980613 




1.096 


0.5 


3.60 x 10 57 




1 


GRB980703 


6891 


0.966 


2.398 


1.26 x 10 s8 




2 


GRB990123 


7343 


1.6 


16.41 


3.0 x 10 s9 


300 


2 


GRB990506 


7549 


1.31 


18.56 


2.08 x 10 59 




2 


GRB990510 


7560 


1.619 


8.16 


1.54 x 10 59 


300 


2 


GRB990712 


7647 


0.434 


11.64 


8.56 x 10 57 




2 


GRB991208 




0.706 


11.2 


2.72 x 10 58 




3 


GRB991216 


7906 


1.02 


67.5 


4.06 x 10 59 


200 


2 


GRB000131 




4.5 


1.5 


3.35 x 10 69 




3 


GRB000301c 




2.03 


1.32 


4.34 x 10 s8 


90 


3 


GRB000418 




1.118 


3.3 


2.5 x 10 58 




3 


GRB000926 




2.066 


10.0 


3.45 x 10 59 


120 


3 



4 *References for beaming factor follow: GRB 990123: Holland et al. 2000; 
GRB 990510: Holland et al. 2000, Harrison et al. 1999; 
GRB 991216: Halpern et al. 2000; 

GRB 000301c: calculated for t b =4.1 days from Bhargavi & Cowsik (2000) 
using Eq 1. of Sari et al. 1999; GRB 000926: Price et al. 2000 

The references for fluxes are as follows: 1. Schaefer(2000); 2. BATSE catalog; 

3. conversion made by applying band-pass corrections as explained in Appendix B. 
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Fig. 1. The results for the Log- normal luminosity function are shown. Right Panel: Results are shown for the full GRB 
sample. The bigger regions in the center enclose the allowed region from the afterglow observations. The region with 
dotted lines corresponds to a run where beaming corrections are applied. The 4 contours on the left side represent the 
region of K-S probability Pk s > 0.01 for the consistency between observed and theoretical number count-flux relation. 
They correspond, from bottom to top (with increasing luminosity) to four models 1-4 respectively of GRB evolution 
described in the text. Left Panel: Results are shown for the GRB sample with GRB duration (in T90) exceeding 2 sec. 
Model 4 of GRB evolution is not shown. 
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Fig. 2. The results for the Schechter luminosity function are shown. The right and left Panels correspond to the full 
sample and the long duration sample (T90 > 2 sec) , respectively. The bigger regions at the top of the figure come from 
the afterglow observations. The region inside the solid (dashed) curves correspond to no (with) beaming correction. 
The smaller regions at the center are from the K-S test for the number count-flux relation. These correspond, with 
increasing photon luminosity, to GRBs evolution models III and IV discussed in the text. The allowed regions for 
models I and II fall below the allowed regions for the models shown. 
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Fig. 3. The results for the scale-free luminosity function are shown. In this model L m j n = 3x 10~ 2 L* and L max = lOOi*. 
The right and left Panels correspond to the full sample and the long duration sample (T90 > 2 sec) , respectively. The 
regions at the left of the figure come from afterglow observations; the region inside the solid (dashed) curves correspond 
to no (with) beaming correction. The thin regions are from the K-S test for the number count-flux relation. These 
correspond, with increasing photon luminosity, to GRBs evolution models III and IV discussed in the text. The allowed 
regions for models I and II fall below the allowed regions for the models shown. 
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Fig. 4. The results for the scale-free luminosity function are shown. In this model L m i n = and i max = lOOOL,. 
The right and left Panels correspond to the full sample and the long duration sample (T90 > 2 sec) , respectively. The 
regions at the left of the figure come from afterglow observations with no (with) beaming correction corresponding to 
solid (dot-dashed) curves. Thin regions in the center of the figure are from the K-S test for the number count-flux 
relation. These correspond, with increasing photon luminosity to GRBs evolution model III and IV discussed in the 
text. The allowed regions for models I and II fall below the allowed regions for the models shown. 



